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Abstract. Macroscopic quantum phenomena (MQP) is a relatively new research venue, with 
exciting ongoing experiments and bright prospects, yet with surprisingly little theoretical 
activity. What makes MQP intellectually stimulating is because it is counterpoised against the 
traditional view that macroscopic means classical. This simplistic and hitherto rarely challenged 
view need be scrutinized anew, perhaps with much of the conventional wisdoms repealed. In 
this series of papers we report on a systematic investigation into some key foundational issues of 
MQP, with the hope of constructing a viable theoretical framework for this new endeavour. The 
three major themes discussed in these three essays are the large N expansion, the correlation 
hierarchy and quantum entanglement for systems of 'large' sizes, with many components or 
degrees of freedom. In this paper we use different theories in a variety of contexts to examine the 
conditions or criteria whereby a macroscopic quantum system may take on classical attributes, 
and, more interestingly, that it keeps some of its quantum features. The theories we consider 
here are, the O(N) quantum mechanical model, semiclassical stochastic gravity and gauge 
/ string theories; the contexts include that of a 'quantum roll' in inflationary cosmology, 
entropy generation in quantum Vlasov equation for plasmas, the leading order and next-to- 
leading order large N behaviour, and hydrodynamic / thermodynamic limits. The criteria 
for classicality in our consideration include the use of uncertainty relations, the correlation 
between classical canonical variables, randomization of quantum phase, environment-induced 
decoherence, decoherent history of hydrodynamic variables, etc. All this exercise is to ask only 
one simple question: Is it really so surprising that quantum features can appear in macroscopic 
objects? By examining different representative systems where detailed theoretical analysis has 
been carried out, we find that there is no a priori good reason why quantum phenomena in 
macroscopic objects cannot exist. 



- This is a later version than the Feb. 24, 2011 version appearing in J. Physics (Conf. Series) for DICE 2010 meeting. 
1. Quantum / classical, micro / macro 

There are many ways to deal with the issue of quantum-classical correspondence [T] . In the most 
common and traditional view the classical limit corresponds to h — > 0, or, invoking the Bohr 
correspondence principle, the principal quantum number of a system n — > oo, or regarding the 
coherent state as the 'most classical' quantum state, or the Wigner function as the 'closest to 
classical' distribution. Less precise criteria also abound, such as the loose concept that a system 
at high temperature behaves classically, or viewing the thermodynamic / hydrodynamic limits 
(of a quantum system) as classical. (For a description of the various criteria, see, e.g., [2]). 



There are holes in almost all of the above common beliefs. A more sophisticated viewpoint 
invokes decoherence, the process whereby a quantum system loses its coherence (measured 
by its quantum phase information) through interaction with its environment [3]. In this 
work we examine an alternative perspective, as the folklore goes, that quantum pertains to 
the small (mass, scale) while classical to the large (size, multiplicity). This common belief 
now requires a much closer scrutiny in the face of new challenges from macroscopic quantum 
phenomena (MQP), viz, quantum features may show up even at macroscopic scales. A common 
example is superconductivity where the Cooper pairs can extend to very large scales compared 
to interatomic distances and Bose-Einstein condensate (BEC) where a large number N of 
atoms occupy the same quantum state, the N-body ground state. Other examples include 
nanoelectromechanical devices [lj where the center of mass of a macroscopic classical object, 
the cantilever, obeys a quantum mechanical equation of motion. Experiments to demonstrate 
the quantum features such as the existence of interference between two macroscopic objects have 
been carried out, e.g., for C 60 molecules passing through two slits [5J or proposed, e.g., for two 
mirrors [6j [7] . 

A most direct account of the difference between the microscopic and the macroscopic 
behaviours of a quantum system is by examining N, the number of physically relevant (e.g. 
for atomic systems, forgetting about the tighter-bound substructures) quantum particles or 
components in a macroscopic object. One may ask: At what number of N will it be suitable 
to describe the object as mesoscopic with qualitatively distinct features from microscopic 
and macroscopic? In the case of classical equilibrium statistical mechanics this issue is 
easily addressed by use of the (grand) canonical ensemble, and such inquiry is answered by 
understanding how the (grand) partition function behaves as a function of N and temperature 
T (and chemical potential fi) all the way to the thermodynamic limit. In classical nonequilibrium 
statistical mechanics, this issue underlies the derivation from molecular Hamiltonian dynamics 
the thermodynamic and kinetic properties (such as transport functions) of a gas of N molecules, 
and their dynamics, which possess salient dissipative and time-asymmetry features nonexistent in 
the microscopic dynamics. Boltzmann answered this question with great valour and magnificent 
success, bringing forth also the issue of time- arrow and providing a molecular dynamics basis 
for the Second Law. To perform a quantitative analysis of such issues one needs to work with 
kinetic and stochastic equations for the N particles, the Boltzmann equation -BBGKY hierarchy 
for (effectively) closed systems in the case where there is no distinguished party, and when there 
is, the Langevin equation for open systems. In the recent decade significant advances have 
been made in providing a molecular dynamics basis to the foundations of thermodynamics [8J, 
relating the macroscopic thermodynamic behaviour of a gas to the chaotic dynamics of its 
molecular constituents. One could even calculate the range in the number of molecules where a 
microscopic system begins to acquire macroscopic behaviour and hence identify the approximate 
boundaries of mesoscopia [§] . 



For quantum systems one needs to deal with additional concerns of quantum coherence and 
entanglement which are critically important issues in quantum information processing (QIP) 
|10| . A fundamental issue in QIP is how the performance of a quantum information processor 
alters as one scales the system up. This dependence on N is known as the "scaling" problem 
We will defer considerations of quantum entanglement of macroscopic objects to a third 
paper in this series after we have a chance to explore how quantum correlations and fluctuations 
impact on MQP in the second paper using the n particle-irreducible (nPI) representation. There 
are many important and interesting issues of MQP, one subset of special interest to us is how 
quantum expresses itself in the macroscopic domain since usually macro conjures classicality. 
Thus even the simplest yet far from naive question need be reconsidered properly. For example, 
why is it that an ostensibly macroscopic object such as a cantilever should follow a quantum 
equation of motion. This center-of-mass axiom is implicitly assumed in many descriptions of 



MQP but rarely justified. The conditions upon which this can be justified are provided by two 
of us and a different third author |12j with the derivation of a master equation for N harmonic 
oscillators (NHO) in a finite temperature harmonic oscillator bath. 

In this paper we ask the question: When will a quantum system of N particles or components 
begin to acquire or show classical features. We shall present an array of physical examples 
to illustrate how this comes about, and glean from them 1) how classicality is defined, 2) 
How the meaning of classicality changes in different types of theories? For 1) the criteria 
mentioned at the beginning and expounded further in the next section will be invoked in several 
concrete examples, such as, how the inflaton dynamics determines when the inflationary universe 
begins to behave classically, from the simplest criteria (invoking the uncertainty principle) to 
the more involved (correlation between canonical variables) and sophisticated (decoherence) . 
Other criteria studied are randomization of quantum phase, entropy generation, decoherence of 
hydrodynamic variables due to the existence of conservation laws. For 2) we give a detailed study 
of a) the O(N) model: an N component quantum scalar field (j) with X(j> self interaction (in 
section [3]) (see, e.g., |13 |. I14| . fT5] ) . the quantum mechanical version of which describes a system of 
N interacting quantum harmonic oscillators, as it is the common underlying structure of many 
physical examples, where some claims of macrocity (being macroscopic, for lack of a better noun) 
and classicality originate. For our purpose here it is sufficient to work with a zero-dimensional 
field theory, i.e., a quantum mechanical model - not only can one shred off the shrouded field 
theory technicalities but it is more directly related to atomic-optical and condensed matter set- 
ups for easier experimental comparisons, b) semiclassical stochastic gravity [16], and c) gauge / 
string theory (in section [5]) . We use some well-known results in quantum field theory - the large 
N expansion - to consider mean field behaviour and under what conditions the semiclassical 
limit of a quantum theory will be reached at the leading order (LO) large N expansion, and the 
role of quantum fluctuations at the next-to-leading-order (NLO) large N expansion. 

Both the O(N) theory and semiclassical gravity show semiclassical behaviour when N becomes 
very large. For the former the leading order large N expansion yields a mean field theory with 
Vlasov equation which is time-reversal invariant. We use this model as an example to expound 
these issues in section [H The next-to-leading order (NLO) expansion begins to show dissipative 
behaviour in the two particle irreducible (2PI) representation [T7]. Using the closed-time-path 
(CTP) 2PI NLO effective action an H-theorem can be shown to exist for the quantum mechanical 
O(N) model [18J. In the latter example, perturbative quantum gravity interacting with an 
N component quantum scalar field emerges as the semiclassical limit when N becomes very 
large [19J. The equations of motion (semiclassical Einstein equation) is time-reversal invariant. 
The NLO expansion yields stochastic gravity [2D] which has both dissipative and fluctuation 
features. These prior results provide a stimulus for us to inquire about the quantum-classical 
correspondence issue in an alternative way, i.e., instead of using the loop expansion (in orders 
of 1/h) where quantum features reveal at successively higher loop orders as corrections to the 
classical solution, we use the large N expansion, where N provides a quantitative measure of 
the 'magnitude' or 'complexity' of a quantum system, to examine how the semiclassical limit of 
a quantum theory is reached. We shall rely heavily on the results obtained in |2H [22l [23] to 
illustrate these crucial aspects in MQP. However, as the examples in gauge and string theory 
show |24[ [25l [26] (section [5]) the leading order large N does not necessarily correspond to a 
classical limit. It may correspond to a thermodynamics limit instead. These are the purposes 
for our invoking different theories to illustrate different criteria and see their fine and not so fine 
distinctions. 

Perhaps it is necessary to make this remark here before we describe the contents of this paper: 
many ingredients in our discussions are well-known to different communities, such as the large N 
expansion in field theory and critical phenomena, dissipation and fluctuations in nonequilibrium 
statistical mechanics, stochastic equations in quantum open systems, hydrodynamic variables 



and quasi-classical domains in decoherent histories, etc. We see this paper's contribution is a 
modest one: discerning the physical essence of apparently similar technology and sharpening the 
often mixed-up physical meanings while focusing on the issues interspersed between the macro 
and the quantum. We hope to use these well-known parts to elucidate some lesser known, even 
misunderstood facts, and to shape a new perspective, as a first step towards constructing a 
theoretical foundation for MQP. [U 

Readers familiar with these theories should just skip over the technical details in those parts 
and concentrate on how they bring out the physical issues raised. In section [2] we display an 
array of criteria behind the common beliefs pertaining to quantum, classical and macro, draw 
the linkages between them, differentiate seemingly equivalent criteria, and issue some warnings 
to over- generalizations. In section Owe present a) the O(N) model, give a quick summary of the 
large N expansion method, leading to a set of coupled equations where one can see the variation 
with N in physical quantities of interest. In section 0] we use three examples to describe how a 
macroscopic quantum system acquires classical features: 1) The so-called 'quantum roll' model 
of inflationary cosmology [27] based on the O(N) model, how stringent the different classicality 
criteria are: amongst them the use of the uncertainty relation, correlation in conjugate variables, 
and environment-induced decoherence. 2) Quantum Vlasov equation (QVE) as an example of 
a mean field theory obtained from the leading order large N (LOLN) expansion of the O(N) 
model. Following [28] we express QVE as a coupled set of equations, one for the number of 
particles, the other for their correlations which contain quantum phase relations. This enables 
one to see how coarse-graining over the phase information brings about entropy generation 
which is related to uncertainty which can be used to quantify classicality. We also mention 
how in the decoherent history conceptual framework quasi-classical domains can appear owing 
to the existence of conservation laws for hydrodynamic variables, thus highlighting the role of 
emergence and importance of collective behaviour in MQP. In section [5] we consider the next- 
to-leading order large N expansion and try to identify the role of quantum fluctuations in MQP, 
using examples given above to make explicit the conditions a macroscopic system can shred 
off or retain its quantum features. We use two more illustrative theories here, b) semiclassical 
stochastic gravity and c) gauge / string theories to examine their LOLN and NLO dynamics 
and explore their classical versus quantum behaviours. We learned that there is no a priori 
reason why a quantum system with large N need to be classical. We discuss the differences 
in physical meaning between a mean field and a classical field, and under what conditions the 
former leads to the latter. Only for the class of Gaussian theories like the O(N) model where 
the mean field theory obtained from the LOLN expansion is equivalent to the classical theory 
would it make sense to think of deviations (from LO expansion) from the mean as providing 
the quantum corrections which then allows one to use large N to address the quantum- classical 
correspondence issue. In section we summarize what we have learned in this intellectual 
exercise and make some general observations. The overall lesson we learned is that we need 
to be very careful in finding the conditions and formulating criteria for a system with large N 
to behave classically. To be safe, one should think of all systems as quantum intrinsically at 
all orders of N, and not to identify mean field as classical without knowing the structure and 
behavior of the theories we are talking about. 

1 Because we want to be as concrete as possible in making our points we decided to extract technical details 
verbatim from some representative papers on the specific topics, instead of grossing over and talking abstractly. 
For this we show our appreciation to the authors of these papers for their lucid presentations which facilitate our 
illustrative purpose. An apology to the readers goes along with our appreciation to the authors: To experts in 
field theory, string theory, gravity and cosmology, the themes in this paper may be too familiar to them and the 
inclusion of details may thus make the contents overladen. This is because we need to consider readers from the 
atomic, optical and condensed matter community where MQP experiments and theories are currently undertaken. 



2. Large N in relation to quantum and classical 

In this program we have set up for ourselves we need to clearly distinguish the relation 
between micro-macro / quantum-classical as well as the thermo-hydro limits. The micro-macro 
correspondence underlies the basic theme of statistical mechanics which has been investigated for 
centuries, and although mesoscopic physics has been with us for three decades many key issues 
bearing on quantum correlations and fluctuations remain a challenge today The quantum- 
classical transition /correspondence has been studied in earnest in the last two decades and we 
have gained much understanding of the role environments play in bringing about a quantum 
system's classical behaviour, but the issue of quantum and macro remains little explored, since 
most common beliefs associate micro with quantum and macro with classical. MQP is a strong 
statement that this need not be the case. We ought to scrutinize all conventional assumptions 
and reexamine every single concept and criterion before reconnecting them into a new conceptual 
framework to meet this new challenge. 

As a background for motivation and a collection of key points we list below the common 
beliefs for a macroscopic object to assume classical behaviour. We pay special attention to the 
relation of large N, quantum and classicality. 

2.1. h, loop expansion and classicality 

• Why is it that % — > is often viewed as classical? This can be easily seen from the path 
integral formalism: 1/fi multiplies the action in the path integral, which has no further 
dependence on h. The dominant contribution to the path integral comes from the path 
which extremizes the action, i.e., the classical trajectory. The steepest descent method 
becomes exact when % — > 0. It is in this sense that taking Ti — > give rise to classical 
physics. 

Suppose the macroscopic object is made up of a space-time inhomogeneous boson 
condensate. Its classical behaviour manifests itself via the vacuum expectation value of 
a scalar field tp: 

<f>{x,t) = (0|V(x,t)|0) 

where |0) is the vacuum with the space-time inhomogeneous boson condensate creating the 
macroscopic object. Note that i\) n contains not only n , but also the product of normal 
product terms. Rearranging a product of normal product terms into a normal product 
creates c- numbers due to the contraction of creation and annihilation operators, while each 
contraction creates a c-number of order h, namely, 

(0|V(x,t) n |0) =<p{x,t) n + 0{K). 

Thus only when the contractions which produce the loops in the Feynman diagrams are 
ignored can one replace (0|V' n |0) by (f> n . This gives the tree level approximation which is 
usually termed "classical" . 

2.2. Large N, mean field and fluctuations 

• Fluctuations to mean: Bohr's correspondence principle states that the classical behaviour 
of a quantum system with n quantum levels is established when the quantum fluctuation 
An is much smaller than the average quantum number n, i.e. An/n <C 1. This condition 
does not depend on h, but is controlled by the number of participating particles or quanta. 

• Large N as a mean field: The limiting theory when N — > oo in a leading order large N 
expansion is usually referred to as a mean field theory in that contributions from fluctuations 
are maximally suppressed. To see this note that in every theory known to have a sensible 



large ./V limit, the vacuum expectation of any product of operators, AB, satisfies the 
factorization relation 

(Ab) = (A)(b) + o(i/n). 

(This means that the disconnected Feynman graphs always dominate). Therefore the 
variance of any operator vanishes as N — > oo, 

ton ((A 2 ) - (A) 2 ) = 0. 

N-¥oo 

This means that the fluctuations become negligible when ./V — > oo. 

• h — > and N — > oo are different limits. The loop expansion and large N expansion are 
very different. In QCD, N determines the number of degrees of freedom (./V 2 — 1 gluons), 
therefore the limit N — > oo resembles the infinite volume limit of a two-dimensional lattice 
system. This is different for theories involving A^-component vector fields such as the O(N) 
model described in section [3j There, the H — > limit gives the complete classical theory 
whereas N — > oo is the gaussian approximation of a quantum theory, for which classical 
and quantum dynamics are identical. So even for this class of theories these two limits are, 
strictly speaking, different. 

2.3. Large N, classical or thermodynamic limits? 

• Is the mean field a classical field? There is a common misunderstanding which equates 
mean field with classical field, whence N — > oo limit is viewed as a condition of classicality. 
But as pointed out by Habib [29, 30j(and reference therein) it is only for Gaussian systems 
such as the O(N) QM model that this is true. It is not true, e.g., for gauge fields, as shown 
below. 

• Large N as thermodynamic limit? Coleman in his 1979 Erice Lectures |31| . perhaps the 
most commonly read introduction to large N field theory, drew a parallel with the classical 
limit: "There is a classical gauge field configuration, which I will call the master field, such 
that the large-N limits of all gauge invariant Green functions are given by their values at 
the master field." Along this line Witten in his 1979 Cargese lectures p5j presented an 
argument for the existence of a master field in the large- N limit of QCD. However, O. Haan 
|32| pointed out that a crucial assumption Witten made, that the expectation value of any 
U(N) symmetric operator of the form f(A) = Ti[A(xi)...A(xi)] where A are N x N matrix 
fields, factorizes, i.e., < f 1 (A)...f k (A) >= IT/ =1 < f t (A) > [1 + 0(1/AT 2 )], does not hold. 
He showed that for the two matrix models this 'masterfield' does not exist. Instead of the 
classical limit Haan argued that the factorization of expectation values of invariant operators 
is a property analogous to the vanishing of fluctuations for macroscopic observables in the 
thermodynamic limit. We will return to this issue in section 5.3. 

• Hydrodynamics and thermodynamics limits. See section 14.41 on this issue. 

2-4- Coherent state and classicality 

• Yaffe in |33| stated that to view Large A" limits as classical mechanics requires 4 assumptions. 
If these assumptions are satisfied, one can generate a natural set of generalized coherent 
states. These coherent states may then be used to construct a classical phase space, derive 
a classical Hamiltonian, and the resulting classical dynamics is equivalent to the limiting 
form of the original quantum dynamics. 

We note that by way of a microscopic model in the open quantum system framework Zurek, 
Habib and Paz [33] studied environment-induced decoherence and showed that there is a 
tendency for a quantum system to evolve to a coherent state which is commonly viewed as 



the 'most classical' quantum state. This is related to the uncertainty relation criterion for 
classicality as illustrated in [21 [35] . 



3. Quantum mechanical O(N) model: large N expansion 

Large N expansion of the O(N) model in 3+1 dimensional quantum field theory has been 
investigated extensively in the literature. The results of immediate relevance for our purpose are 
from the work of Cooper, Dawson, Habib, Kluger, Mihaila, Mottola, Paz et al. For the O(N) 
quantum mechanical model describing the dynamics of a system of N nonlinear oscillators we 
follow the treatment of [2H 122], while relegating the 2PI aspects [181 [23] to our sequel paper. 
For this model numerical integration valid for all N is possible with the quantum roll initial 
condition [27J which provides a measure of how accurate the large N (LO and NLO) expansion 
is as a function of N. We will make use of these results to analyze the quantum-classical problem 
for large N quantum systems. 



3.1. O(N) model 

The Lagrangian for the 0(N) model in quantum mechanics is given by: 

1 N 

L(x,x) = -Y,x?-V(r), (3.1) 
z »=i 



where r 2 = J2i=i x i an d V(x) is a potential of the form 

y(r) = 8?7 ( r2_r ° 2 ) 2 (3 ' 2) 

where g is the self-coupling constant (related to the A in its field theory origin) and vq is where 

the potential goes to zero, V{tq) = 0. The time-dependent Schrodinger equation for this problem 
is given by: 

. dijj(x, t) 



at 



V>(x,t). (3.3) 



The initial conditions for the quantum roll problem allow a numerical solution for all N, as only 
the radial part of the wave function enters. With the exact solution available one can examine 
the validity of the results from a large- N expansion. 
Assuming a solution of the form 

i;(r,t) =r( 1 - N V 2 ( j)(r,t), (3.4) 
the time dependent Schrodinger equation for (j)(r, t) reduces to 



dt \ 2 dr 2 

with an effective one dimensional potential U (r) given by 



* M " f + (3.5) 



(N-l)(N-3) g , 2 2 x2 

u <r) = ^2 + mV ~ r V ■ m 



Upon the rescaling 



r 2 = Ny 2 , r 2 = Ny 2 



(3.7) 



it becomes 



U(y) (N-l)(N-3) g 2 2 2 



N 8N 2 y 2 

corresponding to the new Schrodinger equation, 

f 1 S 2 



dt 



2N 2 dy 2 



+ u(y,N)} <f>(y,t) 



(3.8) 



(3.9) 



where t = N t. 

To implement the large N expansion, following the original work of [37], we introduce a 
composite field \ and add a constraint term 



N 
2^ 



9(2 2\ 

X- ^(r -r„, 



to the original Lagrangian which yields an equivalent Lagrangian, 



L'(x,x,x) = J2o ( j 



n N 2 

■<7 - V'7 ) + yX+ • 



(3.10) 



(3.11) 



One can also glean off the results from the corresponding field theory treatment in |13] (by 



specializing to + 1 dimensions, and replacing (j) a (t) — > Xi(t), /x 2 



r o9 
' 2N 



, and A — > g.) 



3.2. Large N expansion 

The generating function Z[j, J] is given by the path integral over the background fields Xi(t): 
Z[j,J] = e iH ^ J ] = [d X Rdx, e W {iS[x, X ;j,J}} , 

i 

S[x,X',j,J] = J^dt^L' + ^jiXi + Jx) ■ 

The effective action, to order 1/N, is obtained by integrating the path integral for the generating 
functional for the Lagrangian (|3.1ip . over the Xi variables, and approximating the integral over 
X by the method of steepest descent (keeping terms up to order 1/N). A Legendre transform 
of the resulting generating functional then yields the effective action, which we find to be: 

r[g,x] = ^dt|^[gf(t)- X (t)^(t)]+^l n[ G^(i,t)] 



r 2 M i 1 

+fx(t) + Yg x 2 (t) + -HD- 1 (t,t)]j, 



(3.12) 



where the integral is over the close time path C, discussed in [13] and q(t) = (xi(t)). Here 
G7/(i,0 and D~ l (t, t') are the lowest order in 1/N inverse propagators for Xi and Xi given by 



G^(t,t') = \^ + x(t)\WtJ)&ij = Gr x (tJ)8 ij , 



D-\tX) 



N 

-S c (t,t')-U(t,t') 
9 



where 



n(M') = ~E G «(<.OGii(t',t) 

Here 5c(t,t') is the closed time path delta function. 

The equations of motion for the background fields qi(t), to order 1/N, are 



(3.13) 



d^ 
dt 2 



+ X(t) «W + iE/ c dt' Gij (t, t') D(t, t') gj (t') = , (3.14) 



with the gap equation for given by 



(3.15) 



(9) 

The next-to-leading order X{ propag ator Gij(t,f) and self energy Sjj(t,i') to order 1/N are 

4 2) (t,0 = GijM-Y, I d *i / ^Gikit^kiihMG^t') , 

k,l C C 

E kl (t,t') = iG k i(t,t')D(t,t')-q k (t)D(t,t')qi(t') . (3.16) 

(compare these equations with (2.18-2.22) of [13].) 

The actual equation for Q which follows from the effective action differs from (|3.16|) in that 
the final G in the integral equation is replaced by the full Q. This leads to a partial resummation 
of the 1/N corrections which guarantees positivity of (x 2 (t)} (but not the full positivity for the 
density matrix). 

In order to solve for D(t,t'), we first write 



-D(t,t') = -5 c (t,t') + -AD(t,t') 
9 9 



then AD(t,t') satisfies the integral equation, 

N 
9 

(Compare with (2.13-2.16) of Ref. [13].) 



AD(t, t') = |-n(M')- / dt"U(t,t")AD(t",t') , 
N Jc 



(3.17) 



(3.18) 



3.3. Effective potential 

As pointed out by Mihaila et al [21] the effective potential in the large N approximation has 
been previously obtained by Root [381 to order 1/iV. We continue to follow their treatment. 
From the action 13.121 they found that □ 

^(nx)-^^-!)*^ (3-19) 
+ T/£^(*)] + i/£H*-'C*)l, 

2 When Xi and x are independent of time, one can ignore the closed time path ordering and use Fourier transforms, 
passing the poles by using the Feynman contour. 



where x satisfies the requirement 



' i 

7r V eG (r, X ) = 0. (3.20) 
OX 

Using the rescaled y variables defined in (|3,7p the effective potential (|3.19p becomes 

v l$(y>x) x (2 x 2 Vx rs21 s 

1 , ^ x 

+ ("H + m- -3 VX) • 



where m± = 6 db Vb 2 — c, with 



c = ^x 2 + g(^y 2 x + lvx^j (3.23) 

The gap equation which determines x follows from (|3.2U|) 

* = 2 (y - yo) + TiVV^ + 2iV % • (3 - 24) 

To leading order in the large N expansion, (I3T2TD . (13~221) reduce to the parametric set| 

vSix) _ 

TV 2 9 4 ' 

y 2 (x) = yl + ^x-^- (3.25) 

These snugly placed equations are very useful for our purpose of analyzing the aforementioned 
issues in MQP. A result of significance is finding the lower bound on N whereby the NLO results 
give a good approximation to the exact. How does this bear on the issues at hand, namely, 
quantum behavior and classicality in macroscopic systems? We will discuss this in section 5.1. 

In the next section we will focus on the leading order results which are much simpler but still 
provides useful insight on the macro-quantum relation. 



4. Number - phase / correlation, coarse-graining / uncertainty - entropy, 
decoherence / classicality and hydrodynamics 

In this section while staying at the leading order large N (LOLN) level we adopt another angle 
to analyze mean field theory and classicality. Using the uncertainty relation to find out when 
a quantum system begins to show classical behaviour is shown to be inadequate, as we will 
illustrate with an example of quantum roll in the inflationary universe. Decoherence and 
correlation are needed. In the second example we use a theory shown to be obtainable as 
the LOLN of an O(N) model which is time-reversal invariant and examine another set of criteria 
for classicality, that of quantum phase information and entropy generation. Although a system 
with no quantum phase information is usually identified as classical, yet in essence this is not the 

3 These equations (|3.21[) and p.24|) agree with Root [38], however he used the leading order expression for \ i n 
IpPgJl . rather than the full \ of (EH}. 



case. Fluctuations in the number density carry phase information which signifies its quantum 
origin. One needs to ensure that phase information at all orders vanish for all times to be able to 
say so, but this is not natural and not easily implementable. In the third subsection we discuss 
classicality from the decoherent history viewpoint, that hydrodynamic variables are more readily 
decohered because of the existence of conservation laws for these variables. We also mention 
why thermodynamic and classical limits are not identical. 

4-1- Criteria of classicality: uncertainty relation, decoherence and correlation in a quantum roll 
One commonly used signifier of when a quantum system begins to show classical features is the 
use of the uncertainty relation to distinguish these two regimes. (See [2].) A concrete example 
with some importance in early universe cosmology is Guth and Pi (GP)'s [27] use of a quantum 
mechanical model for the description of the inflaton dynamics of a scalar field (the inflaton). 
They solve for the evolution of the wave function in a de Sitter universe to illustrate how 
the inflationary transition takes place. The O(N) model has been used in Friedmann universe 
models by Boyanovsky [39], Cooper et al [30] and many others. The scalar field starts at the 
top of an inverted harmonic oscillator potential and 'rolls down'; this so-called 'quantum roll' 
initial condition is in fact used in the results reported in the previous section by Mihaila et al. 
We want to use these studies to consider the relation between large N, macroscopic quantum 
systems, the uncertainty relation and classicality. We first give a brief account of how GP 
treated this problem, then point out the deficiencies in the criterion they used for classicality. 
We then describe the necessary requirements for classicality as represented by the behaviour of 
the Wigner function, namely that it decoheres in the presence of an environment and it possesses 
correlations between the canonical variables as in a classical trajectory in phase space. 
GP considered an eternal inflation situation with the background metric in the form 

ds 2 = -dt 2 + a(t) 2 dx 2 (4.1) 

where the scale factor a(t) = e Ht for eternal inflation with the Hubble expansion rate given by 
H{t) = ^/87rGpo/3 where po is the vacuum energy density of the false vacuum. The inflaton 
field (j) dynamics is described by the action 

S = J d A x^-g{- IgTd^t - V {$)] (4.2) 

with the potential V((p) given by 

V^) = \\l^ 2 -^ (4.3) 

where \x is related to the value of the field (j) c = ^= at the minimum of the potential. The false 

vacuum energy is given by po = j\, giving an expansion rate H = ( 2 ^-)zp 2 . From the action 
the inflaton field obeys the dynamical equation 

4> + 3H<j) - a - 2 (t)A^(j) = -dV/d(f) (4.4) 

where A^ 3 ) denotes the 3-dim Laplace operator on the spatial hypersurface. 

GP then expanded the free field (assumed to be confined within a 3-dim box of length L) 
in Fourier modes with o~ K being the amplitude for the ± propagating components of the A;th 
mode and the zero mode (treated separately because of infrared divergence considerations) and 
its conjugate momentum ir K . 



As initial condition GP assumed that at early times the system was in thermal equilibrium 
with a heat bath at temperature To which rapidly drops as the universe inflates away T = Toe~ Ht . 
The expectation value of physical quantities such as a K with respect to this thermal ensemble 
can be calculated by taking the ensemble average, 



where ip is given by the Hankel functions of the first kind and Q K is the initial value of j^ry at 
t — > -co. 

To verify that the system behaves classically at late times GP calculated how the uncertainty 
function changes in time. 



The denominator h/2 is the value for minimal uncertainty in the case of a Gaussian wave packet 
at zero temperature. Physically the k mode amplitude can be viewed as behaving classically 
when its effective wavelength \ e jj = 2ir/k e ff (where k e ff = a(t)/p(t), with a(t) = e Ht , 
p K = \J\& R + 7 2 , 7 = \J \/TiTq/2 ) is much larger than the horizon length. 

Hu and Zhang [2] (see also Anastopoulos and Halliwell |35j) have derived an expression for the 
uncertainty relation at finite temperature using the Wigner function of an oscillator interacting 
with a thermal bath modelled by N harmonic oscillators as weighing function. They also 
discussed how the contributions of quantum fluctuations weighted against thermal fluctuations 
and the criterion of classicality and thermality. 

As later investigations showed this criterion used by GP for classicality is suggestive but 
incomplete. Two conditions need be satisfied by the (reduced) Wigner function from the 
equivalent (reduced) density matrix ('reduced' refers to the situation where the system is coupled 
to an environment and rendered open). Wigner function has long been viewed as the quantum 
object 'closest to' the classical distribution function. Up to the early 90's before decoherence 
was considered in earnest it was viewed that if the Wigner function of a quantum system peaks 
along a classical trajectory, meaning that there is correlation between the coordinate x and 
momentum p, then one could view the system as classical. This criterion is shown to be valid 
only for Gaussian systems (such as a free harmonic oscillator with quadratic potential coupled 
bi-linearly to an environment of simple harmonic oscillators) because otherwise it may not even 
be positive definite |41j . In the Gaussian approximation the Wigner function is positive for 
all times but as was further pointed out [42J it does not describe classical correlations unless 
the system is coupled to an environment. Thus classical correlation and environment-induced 
decoherence are two necessary conditions for a quantum system to show classical behaviour. 

When the wave packet spreads out further away from the top of the potential the Gaussian 
approximation breaks down. As shown by Lombardo et al (LMM) [43] . as the coupling 
between the system and the environment increases, its decoherence time decreases. Due to the 
nonlinearities of the potential, when the coupling vanishes there is no classical limit, not even 
classical correlations. In fact, as reported by Antunes et al [H], in a quenched phase transition, 
even after classicalization has been reached, the system may display quantum behaviour again. 

With these new understandings LMM reanalyzed GP's model with and without an 
environment, using the quantum Brownian motion |4"6] model, especially results from the 
non-Markovian master equation for a general environment |47| . They found that when the 
system is isolated, due to the high squeezing of the initial wave packet x and p become classically 
correlated. The density matrix is not diagonal. The correlation time depends on the shape of 
the potential. Only when the particle is coupled to an environment will a bona fide quantum to 
classical transition occur. The Wigner function becomes peaked around a classical trajectory 




(4.5) 





(4.6) 



h/2 



and the density matrix diagonalizes. The decoherence time depends on the diffusion coefficient 
in the master equation. 

Now that we have analyzed the quantum mechanical model in detail how would a macroscopic 
system as modelled by the O(N) quantum field theory behave? Symmetry breaking in a 
<fi A 0(N) model has been studied in detail by many groups |14^ I39j. Lombardo et al [13] carried 
out numerical analysis of this model and made the following observations: at late times the 
long- wavelength modes k\ -C t^ 1 in the wave functional become classically correlated. This 
is analogous to the situation described above for the inverted harmonic oscillator. For these 
modes the width of the Gaussian wave function increases linearly with time. Moreover, as the 
width of the wave function or the density matrix increases, the Wigner function becomes sharply 
peaked around the classical trajectory. We see that in the large N limit, and at long times, the 
dynamical evolution of the O(N) model shows classical correlations but not true quantum to 
classical transition. Just as in the inverted oscillator without environment, the density matrix 
does not become diagonal. The correlation time depends on the details of the potential. 

In conclusion Lombardo et al |43j observe that in order to get a Wigner function that is 
positive definite and peaked around a classical trajectory at long times, it is necessary to have 
both vanishing effective mass and a Gaussian wave function. Therefore, it is quite possible that 
in a field theory calculation with finite N there will be no classical limit nor classical correlations 
unless one allows the field to interact with an environment. On this case for an analysis of 
how the nonequilibrium evolution of a condensate and its fluctuations vary with N, Baacke 
and Michalski [48] have provided analytic expressions for general N and numerical solutions for 
N=l,4, 10. This should be compared with the analytic expressions for the effective potential 
calculated up to NLO for finite N given by Mihaila et al for an O(N) model, which we will 
discuss in section 15.11 

4-2. Criteria of classicality: coarse- graining of quantum phase, entropy generation 
We next use the entropy function derived from the number of particles and their quantum 
correlation in a quantum Vlasov equation to illustrate the notion of classicality. The large N 
aspect enters because quantum Vlasov equation is obtained as a mean field theory derivable 
as the leading order large N approximation from scalar QED, where N is the number of 
identical copies of the charged matter field placed in an electromagnetic (EM) field . This 
is sometimes referred to as the semiclassical limit since the matter field is fully quantized and 
the EM field is treated classically. Kluger, Mottola and Eisenberg [28] considered the case of a 
spatially homogeneous electric field, represented by the vector potential in the Coulomb gauge, 
A = A(t)z, Aq = , whence the electric field is given by E = —Az = Ez . 

Assuming also that the field lives in a finite large volume V we can expand the charged 
scalar field operator in Fock space in Fourier modes. Since particles are physically distinct from 
antiparticles, we need two independent sets of destruction operators 



Denote the time-independent annihilation operator of a particle in mode k by and the 
creation of an anti-particle in mode — k by 6^. They obey the commutation relations 




(4.7) 



(4.8) 




(4.9) 



Therefore N + (k) = (a^) , AL (k) = (b^bk) are the mean numbers of particles and antiparticles 
respectively. Without loss of generality we can make use of the freedom in defining the initial 



phases of the mode functions to set the correlation densities (a k a k ) = (^k^k) = 0. In a 
Hamiltonian description we can take for each mode k 



Mt)^h(tW + mt)bl k (4.10) 
as the (complex) generalized coordinates of the field $ and 

Mt) = <&(*) = /k(*)4 + A(*)&-k , (4.11) 

as the momentum canonically conjugate to it. By virtue of the commutation relation (|4,9p they 
obey the canonical commutation relation, 

bk,vr k /] = iM kk / , (4.12) 

provided that the mode functions satisfy the Wronskian condition 

(fkj* k ) = ih. (4.13) 

The complex amplitude function fk(t) of the k th mode satisfies the wave equation 

^ + ul(t)f k (t) = 0, (4.14) 

where the time dependent frequency is given by 

w£(t) = (k - eA) 2 + m 2 = (k z - eA{t)f + k\+m 2 . (4.15) 

where k z is the constant canonical momentum in the z direction while the physical (gauge- 
invariant) kinetic momentum is given by 

Pz(t) = k t -eA(t), p z = -eA = eE (4.16) 

For a spatially homogeneous electric field (i.e., V-E = 0), by Gauss' Law, the mean charge 
density must vanish, i.e, j°(t) = 0. and the mean current in the z direction is 

j(t) = 2eJ d 3 k [k z - eA(t)]\Mt)\ 2 (l + N + (k) + iV-(-k)) (4.17) 

One can further restrict to the subspace of states for which 

N + (k) = AL(-k) = iV k (4.18) 

Clearly the vacuum A^ + (k) = iV_(— k) = (as well as a thermal state) belongs to this class of 
states. 

Particle pairs will be produced in a strong background field, and in turn, affect the strength 
and evolution of this background field. This first step is called a 'test field' approximation 
where the background field is assumed fixed. The second step of including the effects of created 
particles is called a 'backreaction problem'. It demands a self-consistent solution of the mean 
electric field E(i) coupled to the expectation value of the current j(t) of the quantum charged 
scalar field y k (t). In a spatially homogeneous electric field, the only nontrivial Maxwell equation 
is simply, 

- E(t) = A(t) = j(t) (4.19) 

where the current is given by (|4.17p . Since the charged scalar field depends on the vector 
potential A to begin with, / k (t) and A(t) need be solved self-consistently from (|4.14p with 
flUgD and (ITOIl . 



4-2.1. Particle production In a time dependent background the equation for the mode functions 
(|4.14p would in general admit time dependent solutions. This means that one set of particle 
states which define a Fock space at one moment of time would change at a later time into 
another set of particle state in a different Fock space. In particular, an initial vacuum could 
evolve into an n-particle state later, which means particles are being produced. This is best 
expressed in terms of the Bogoliubov transformation [49j between these two Fock space basis ak 
and ak. Consider for the moment just particles, no anti-particles, and no particle interactions 
(see [SD1 EI] for that case). Let the first basis be associated with modes (/k, /£) , the second 

basis ak with modes ( /k, f^j ■ We may expand the field operators in either base, such as 

(t) = h (*) ak + /k (*) al k (4-20) 

in the first case, and 

(t) = h (t) a k + /£ (t) a f _ k (4.21) 
in the second. Since both sets of solutions of the mode equations are complete, we must have 

/k (t) = a k / k (t) + Afc/k (*) , (4.22) 

and its inverse 

/k(t) = «k/k(t) - /3k/k(0- (4-23) 

The Wronskian condition (/k, /k) = (/k> /k) = ^ imposes a condition on the Bogoliubov 
coefficients 

|«k| 2 - |Ac| 2 = 1 (4.24) 

for each k. We can thus write 



|a k (i)| = coshr k (t), 

|/3 k (t)| = sinhr k (t). (4.25) 

where r k (i) is called the squeeze parameter for mode k a terminology adopted from quantum 
optics. For a description of particle creation in the squeezed state language, see, e.g., Chapter 
4 of [52]. 

The above only describes how one set of Fock states defined at one time is related to another 
set, but whether a vacuum can indeed be defined at any one time is a different and often more 
challenging problem. We won't go into the details of this important issue, but just be satisfied 
for our present purpose with the fact that when the background varies sufficiently slowly one 
can invoke the concept of an adiabatic number state 0. This was used by |28j to calculate the 
charged particle production in an electric field problem described above. 

The adiabatic number state /i ) (f) corresponds to the 0th order adiabatic vacuum [defined 
in (4.46) of ED as /+„, ] 

4 This definition of a number state makes use of the fact that under adiabatic evolution, particle number is an 
adiabatic invariant, thus restricting its validity from the start to extremely weak or slowing varying background 
fields. This level of approximation will not give a good measure for on-going particle creation, as particle 
creation is basically a nonadiabatic process. It is however useful for quantum kinetic theory descriptions, where 
a quasi-particle approximation is usually introduced which amounts to incorporating only the quantum radiative 
corrections to the particles but not the fully field theoretical effects such as particle creation. 



where 0^ = J* wi c (t')dt / , is the Oth order adiabatic phase. At this level of accuracy one 
measures particle numbers at all times with respect to the initial vacuum state at time to- The 
adiabatic particle number is defined to be [28] 

iV k (t) = (4(*)Sk(t)> = (&t k (t)L k (t)) = |a k | 2 (4a k > + |/3 k | 2 (6_ k 6L k ) 
= (l + |/3 k | 2 )iV + (k) + |/3 k | 2 (l + iV_(-k)) 

= |/? k | 2 + (1 + 2|/3 k | 2 )iV k = iV k + (1 + 2iV k ) |/3 k (t)| 2 (4.27) 

where the last line is valid only if the number of positive and negative charges are equal (cfr. 
(|4.18p ). The amount of particle production at time t in this basis is given by the expectation 
value of the number operator a^a at time t with respect to the vacuum state | >o defined at to, 
(not the vacuum state | >t defined at t). 



4-2.2. Number and correlation In this time-dependent particle number basis an equation can 
be obtained for the time rate of change of the number of particles created in each mode by 
differentiating (|4,27p . 

^iV k = 2(l + 2iV k )Re(/3 k /3 k ). (4.28) 

We need an expression for /3 k in terms of a k ,/? k and @ k °\t) = J*w k (t')dt / (henceforth we 
will omit the superscript on 6 k °^)- They are given by 

«k = |^/3 k exp(2i9 k ), /3 k = ^-a k exp(-2i6 k ). (4.29) 

and thus, 

iLjy k = ^(l + 2AT k )Re{a k ^exp(-2ie k )} 
at w k 

= — Re{C k exp(-2iG k )} , (4.30) 

where we have defined the time-dependent pair correlation function, 

Ck(i) = (ok(t)6- k (t)) = (1 + 2iV k ) a k /3 k . (4.31) 

The pair correlation C k (t) is a very rapidly varying function, since the time dependent phases 
on the right side of (|4.3ip add rather than cancel. The phases, however, nearly cancel in the 
final combination of (|4.30p to render iV k a slowly varying function. The time derivative of the 
pair correlation function is given by, 

d7 Ck = lir (1 + 2Nk) ex p( 2i0k ) i 1 + 2 i^i 2 ) 



2uj k 



l + 2iV k exp(2iG k ). (4.32) 



4-2.3. Quantum Vlasov equation Let us return now to the two equations for the rates of change 
of the particle number and the quantum correlations. Solving (|4.32p formally for C k , assuming 
that C k vanishes at some t = to which could be taken to — oo, and substituting into (|4,30p we 
obtain 

d^ k = H ft dt ' (* + 2i ^ k(t/) ) cos [20k(t) " 20k(t ' )] } ' (4,33) 



()4.33p gives the rate of particle creation in an arbitrary time-varying mean field which may 
be called a "quantum Vlasov equation". Note the appearance of the Bose enhancement factor 
(1 + 2iVk) in (|4.33p indicates that both spontaneous and induced particle creation are present. 
One important feature of (|4,33p is that it is nonlocal in time, the particle creation rate depending 
on the entire previous history of the system. Thus particle creation in general is a non-Markovian 
process [531 154"! l55| [56] . Note that the nonlocal form of (|4.33p results from solving one variable 
C in terms of the other N, each obeying a Hamiltonian equation of motion. This is a general 
feature of coupled subsystems. 

(|4.33|) becomes exact in the limit in which the electric field can be treated classically, i.e 
the limit in which real and virtual photon emission is neglected, and there is no scattering. 
Inclusion of scattering processes leads to collision terms on the right side of (|4.33p which are 
also nonlocal in general. This nonlocality is essential to the quantum description in which phase 
information is retained for all times. The phase oscillations in the cosine term are a result of the 
quantum coherence between the created pairs, which must be present in principle in any unitary 
evolution. However, precisely because these phase oscillations are so rapid it is clear that the 
integral in (|4.33p receives most of its contribution from t' close to t, which suggests that some 
local approximation to the integral should be possible, provided that we are not interested in 
resolving the short time structure or measuring the phase coherence effects. The time scale for 
these quantum phase coherence effects to wash out is the time scale of several oscillations of 
the phase factor @k(t) — €>k(0> which is of order r qu = 2^/^ = lirh/e^, where ek is the single 
particle energy. 



4-3. Density matrix and entropy 

After the elimination of the rapid variables Ck defined in (|4.3ip in favour of the slow variables Ak 
one can construct the density matrix in the adiabatic number basis easily [25]. In a pure state 
the only nonvanishing matrix elements of p are in uncharged pair states with equal numbers of 
positive and negative charges, 4 = = , with 4 the number of pairs in the mode k, viz. 

<24|p|2^ k ) pure = eKV-'0*k(*) sec h 2 r k (i) (t a nhr k (t))^ k (4.34) 

where the magnitude of the Bogoliubov transformation, r k (t) is defined in (|4.25p and its phase, 
#k(f) is determined by 

a k /3£e~ 2i0k = -sinhr k coshr k e i,?k . (4.35) 

Hence the off-diagonal matrix elements £' ^ £ of p are rapidly varying on the time scale Tq U 
of the quantum mode functions, while the diagonal matrix elements £' = £ depend only on the 
adiabatic invariant average particle number via 



(24 1 p| 2£ k ) purc = p 2 e k = sedrV k tanh 2 Vk 

|/3k| 2 ' k A/j 



k 



k 



(1 + |/?k| 2 ) £k+1 (1 +Ak) lk+1 



(4.36) 



and are therefore much more slowly varying functions of time. The average number of positively 
charged particles (or negatively charged antiparticles) in this basis is given by 

oo 

E 4P2* k = A/k • (4.37) 

£ k =0 

Thus the diagonal and off-diagonal elements of the density matrix in the adiabatic particle 
number basis stand in precisely the same relationship to each other and contain the same 
information as the particle number Ak and pair correlation Ck respectively. 



Prom the (admittedly somewhat simplistic) perspective of the diagonization of the (reduced) 
density matrix as a signifier of classicality, this is similar to the criterion used in environment- 
induced decoherence alluded to in section 4.1, with one major difference: Here, quantum phase 
in the system is eliminated by choice of representation, while there, it is through the system's 
interaction with the environment. The former case needs physical justification for making such 
a choice, the latter needs explicit demonstration of how the pointer basis from measurement 
readings depend on the interaction Hamiltonian. We shall continue the first route below. 

4-3.1. Classicality from ignoring the quantum phase In the density matrix (|4.36p the diagonal 
elements p2i k may be interpreted (for a pure state) as the independent probabilities of creating 
£k pairs of charged particles with canonical momentum k from the vacuum. This corresponds to 
disregarding the intricate quantum phase correlations between the created pairs in the unitary 
Hamiltonian evolution. When physics is expressed in the adiabatic particle number basis (the 
Fock or N representation) the phase information is ignored. The quantum density matrix in this 
representation produces an entropy function which reflects the entropy associated with particle 
creation alone but says nothing about the evolution of the quantum phase or correlation. This 
illustrates the crucial role played by the choice of representations in the definition of entropy 
associated with particle creation [57J. 

Results obtained from neglecting quantum phase are known to be quite accurate for long 
intervals of time in the back-reaction of the current on the electric field producing the pairs, 
because when the current is summed over all the k modes, the phase information in the pair 
correlations cancels very efficiently. Thus for practical purposes one can approximate the full 
Gaussian density matrix over large time intervals by its diagonal elements only, in this basis. 

4-3.2. Entropy generation from particle creation Let us examine the reduced von Neumann 
entropy constructed from the diagonal density matrix (|4,36p 

oo 

SV(*) = -£E^ k hip24 (4.38) 
k e k =o 

Upon substituting (|4.36p into this, the sums over are geometric series which are easily 
performed. The von Neumann entropy of this reduced density matrix 

SAf(t) = £{(1 +A/k)hi(l +A/" k ) -7V k ln7V k } (4.39) 

k 

is precisely equal to the Boltzmann entropy of the single particle distribution function A/"k(t). 
Hence 

^X>(^)^ 

increases if the mean particle number increases. This is always the case on average for bosons 
if one starts with vacuum initial conditions, since |/?k| 2 is necessarily nonnegative and can only 
increase if it is zero initially [58J. Locally, or once particles are present in the initial state, 
particle number or the entropy (|4.4U|) does not necessarily increase monotonically in time. 

Hence the notion of entropy associated with particle creation, and the lore that it increases 
in time, is only valid for spontaneous production of bosons from an initial vacuum state. This 
function associated with fermions, and that associated with stimulated production of both boson 
and fermions, can decrease in time. 



4-4- Decoherence of hydrodynamic variables, conservation laws and quasi-classical domains 
To end this section on the criteria of classicality we want to describe an important formulation 
in how classicality emerges. This is the decoherent histories approach of Gell-Mann and Hartle 
|59| applied to the quantum mechanics of closed systems, which is very different from the 
above-mentioned environment-induced decoherence scheme applied to open quantum systems. 
A central theme of this program is the decoherence of hydrodynamic variables, namely, that 
the variables typically characterizing the quasiclassical domain of a large and possibly complex 
quantum system that will become classical 'habitually' are the integrals over small volumes 
of locally conserved densities. Examples of local densities are number density, momentum 
density, energy density, charge density, which are collectively called hydrodynamic variables. 
The question of emergent classicality then consists first of understanding why these variables 
enjoy this special status, and second, of deriving the familiar hydrodynamic equations for these 
densities from decoherent histories. The representative references are |60t [6T| [62], 163] . The 
treatment in Brun and Halliwell [63J contains examples for ease of understanding. For deeper 
conceptual issues related to quasi-classical domains read Dowker and Kent |64j . 




We mention here those results which are directly relevant to our discussion of hydrodynamics, 
thermodynamics and classicality. 

1) Following Gell-Mann and Hartle [60J Hartle, Laflamme and Marolf [61] showed that exactly 
conserved quantities are exactly decoherent, and it is conceivable that the probabilities for such 
histories would peak about deterministic evolution equations. Brun and Halliwell [63] derived 
a formula which shows the explicit connection between local conservation and approximate 
decoherence. 

2) Halliwell [62] showed that a) for the case of the diffusion of the number density of a 
dilute concentration of foreign particles in a fluid, for certain physically reasonable initial states, 
the probabilities for the histories of the number density are strongly peaked about evolution 
according to the diffusion equation, b) When the initial state is a local equilibrium state, 
large-N statistics guarantees that the probabilities for histories of local densities are strongly 
peaked about a single history which describes the evolution of the mean in the local equilibrium 
initial state. This intuitively clear result readily connects with previous work, which shows that 
the mean values of local densities evolve according to hydrodynamic equations. The fact that 
there is essentially only one history with nonzero probability means that there is approximate 
decoherence. 

3) Brun and Halliwell [63] considered a class of models consisting of a large number of 
weakly interacting components, in which the projections onto local densities may be decomposed 
into projections onto one of two alternatives of the individual components. One model they 
gave is the placement of N particles into two partitions, the other model consists of a long 
chain of locally coupled spins for which the Hamiltonian conserves the total spin. They 
computed the decoherence functional for histories of local densities, in the limit when the 
number of components is very large, they find that decoherence requires two conditions: a) 
the smearing volumes must be sufficiently large to ensure approximate conservation, and b) the 
local densities must be partitioned into sufficiently large ranges to ensure protection against 
quantum fluctuations. 

The last two points are of special interest to us: namely, 2b) referring to the evolution of 
the mean in the local equilibrium initial state shows the relation of thermality and mean field 
dynamics. 3a) involving a large number of components N in a system and a large volume 
V smacks of thermodynamical limit, except that hydrodynamics ensures that conservation laws 
apply for these collective variables. Point 3b) above gives the condition for quantum fluctuations 
to be suppressed. We saw similar descriptions for the mean field and fluctuations in terms of 
leading order and NLO expansions for a large N quantum mechanical system. 




5. Next-to-leading order large N, quantum fluctuations 

So far we have discussed the leading order theory of an 0(N) model obtained from a large 
N expansion. This mean field theory and its behaviour is usually referred to as "classical". 
Assuming such a connection which we will examine critically in this section, we may ask: 1) 
Can one identify a number A c whereby a system containing iV components with A > A c will 
begin to acquire classical features? 2) Since a classical system is devoid of quantum fluctuations, 
can one extract the information related to quantum fluctuations from higher order expansions 
in large N? Following this logic one may hope to come up with a systematic way to identify the 
demarkation zone between classical and quantum features in terms of N for a quantum system 
with a large collection of components. 

But as we forewarned in the Introduction, this equivalence of mean field with classical features 
is valid only for Gaussian theories. More accurately one should think about mean field rather 
than classical theory in the leading order large N and deviations from the mean in the NLO 
large N. One should therefore ask, 3) Does one see a change of qualitative behaviour in the NLO 
or higher order large N expansion from the mean field theory? What measures deviation from 
mean field? Perhaps one should not think of quantum fluctuations, but correlations, as in the 
BBGKY hierarchy? (The role of 2PIEA in closed (democratic) systems such as systems with 
N components will be discussed in our next paper.) Finally, to see the difference of behaviour 
from the O(N) model we have studied above one should also ask, 4) Would the large N limit 
of a gauge field theory behave differently, and how so? What can one say about deviation from 
the mean field and how the contribution of quantum fluctuations enter in these theories? 

In the following three subsections we will address issues raised above: subsection 15.11 in the 
context of O(N) theory for question 1), subsection 15.21 in the context of stochastic gravity for 
queries 2) and 3), and subsection 15.31 in the context of large N gauge field theory for question 
4). 

5.1. How large an N will NLO give sensible results? Quantum features beyond mean field 
theory ? 

For this question let us return to the results reported in section [3] pertaining to next-to-leading 
order (NLO) large N expansion in an O(N) model. The authors of |21J show that the naive large 
A expansion violates unitarity (or more generally, positivity) leading to an instability at least 
for A less than some value Nt- They have provided numerical evidence that a sharp threshold 
at A ~ Nt exists, where for N > Nt they saw no instability. They trace this behaviour to 
be related to the nature of the effective potential at NLO: At this order, the effective potential 
has the property of not being defined everywhere for values of A < N c , where A c depends on 
the values of the parameters specifying the potential, and Nt ~ N c . However, for A > A c , the 
effective potential is defined everywhere and there is no instability. 

Probing further using the effective potential for the NLO expansion these authors asked the 
following question: "How does the point y m in, below which the next-to-leading order large A 
effective potential does not exist, changes as a function of A?" We know that at "infinite A," 
(leading order), y m i n = 0, but it is important to know how this limit is reached. For instance, 
is there a finite value of A beyond which y m ; n = 0? They plotted (figure 5 of [21]) y m ; n as a 
function of A for the set of parameters chosen in their equation (7.15): At the next-to-leading 
order large A level, for A < 18.6, y m i Q is finite, but for A > 18.6, it hits the origin. Thus 
for A > 18.6, the authors concluded, one can associate a quantum state (though not known 
explicitly) with the next-to-leading order approximation. 

These findings by [21] are interesting by implications. If one adheres to the conventional 
wisdom based on the loose concept that mean field theory obtained from LO approximation 
describes classical physics then one might be tempted to think of the NLO expansion as 
containing some quantum features. Does A > A c signals a transition from quantum to classical? 



5.1.1. Is there a critical N c which marks the appearance of MQP? If N = oo is used as a 
criterion to view an N component quantum system as reaching its classical limit, it seems logical 
to ask, "What is the smallest number N c for such classical behaviour in this quantum system 
to disappear? The related question for our concern is, "Does the number N c obtained by these 
authors carry such a meaning, as a demarkation between classical and quantum behaviors?" If 
so this is a very neat result because it provides a quantitative measure in how 'large' N needs to 
be to mark the appearance of MQP. We shall explain in what follows that unfortunately this is 
not the case. Life would have been much easier if it were. The key observation is that it is only 
for Gaussian theories that the mean field theory from LO LN approximation is equivalent to a 
classical theory. The critical N c is only a measure of the validity of the NLO approximation, not 
a measure of the importance of quantum contributions. Quantum features are pervasive at all 
levels of approximation. We will explore this issue further in the context of critical phenomena 
in the sequel paper. 

Convexity and positivity. Let us look at some technical features of the effective potential 
method which these authors used. The generating functional of connected Greens functions 
[W ~ In Z) is convex. This follows from the positivity of the two point function which is obtained 
by functionally differentiating W twice. Legendre transformation takes a convex function to a 
concave function. Paying attention to the signs in the definitions it follows that the effective 
potential is also convex. As a result the convexity of the effective potential is a consequence (or 
a guarantee, depending on which way one looks at it) of the positivity. This is the reason why 
the critical values Nt and N c are about the same. Recall that Nt is the critical value of N 
below which an instability occurs (see figure 7 in MACDH for a plot of < r 2 > as a function of 
t) at finite time and N c is the critical value of N below which the effective potential ceases to 
exist everywhere. 

Perturbative expansions of the effective potential using the background field method run into 
difficulties for classical potentials (tree level) that are non-convex, which is the case described 
in MACDH. The loop expansion develops an imaginary part because of the unstable modes, 
which prompt tunneling or nucleation. If one were to sum up all the loop corrections the result 
would be a convex effective potential (with minima at the origin) . However at finite order there 
usually are domains where the effective potential becomes imaginary. 

The NLO LN expansion is not immune to this shortcoming. It sums up graphs from all orders 
in loop expansion but still the dynamics is altered. The same problem as in loop expansion 
occurs here as well. For N > N c convexity is restored. The effective potential is convex for both 
classical and quantum O(N) models if we don't make any approximation. Hence the fact that 
convexity is restored for N > N c is a statement about the validity of the 1/N approximation, 
not about the system becoming classical or if enough quantum details have been retained. 

5.1.2. Large N approximation does not capture full quantum dynamics As MACDH mention in 
their paper although NLO is a significant improvement over LO for sufficiently large N, at late 
times NLO fails to capture the nonlinear effects as well as LO does. Hence having N > N c (or 
N > Nt, since they are both of the same order) doesn't mean NLO approximation describes 
the dynamics accurately, it simply means the divergences and nonpositivity is avoided. At late 
enough times the approximation will fail (see figure 1 1 in MACDH) . In general to make a claim 
about the validity of an approximation time need be specified, since any approximation will fail 
after some time. To inquire about the turnover from quantum dynamics to classical dynamics 
we are referring to the full dynamics, not just for short time transient or long time asymptote 
behaviours. 

Thus, we see that N > N c is & criterion for the validity of the NLO approximation, as these 
authors meant it to be, not so much about the quantum versus classical behaviour of the N 



component system. As a side note: it is the NLO approximation that has this pathology not the 
LO. LO approximation most likely gets inaccurate and unreliable before the NLO approximation 
diverges. So LO approximation behaving well doesn't say anything about its accuracy. Same 
applies to NLO approximation for N > N c . 

5.2. Semiclassical and stochastic gravity as LO and NLO large N theories 

We want to use these two theories as an instructional example to two important themes in our 
program, 1) the autocratic open quantum system (Langevin) paradigm, where gravity is treated 
preferentially as the system of special interest over the matter field regarded as its environment 
where only its coarse-grained (or averaged or integrated over) effects are accounted for. This 
is how the Einstein-Langevin equation (ELE) [65] was first derived. 2) the democratic large N 
(Boltzmann) system. We shall use an N component matter field as source to illustrate how the 
semiclassical Einstein (SCE) equation and the ELE can be derived as the leading order O (1) 
and next to leading order 0(1/N) theories in the large N expansion respectively. 

Semiclassical gravity refers to the theory where gravity (metric function g ab ) is treated 
classically but matter as a quantum field (f>j, the vacuum expectation value (vev) of whose 
(renormalized) stress energy tensor serves as the source to the semiclassical Einstein equation 

G ab [g] = K (f ab [g}) . (5.1) 

\ / ren 

where the subscript ren denotes renormalized. 

Stochastic gravity [16] includes in the source of the SCE equation an additional term £ ab 
accounting for the fluctuations of the stress tensor operator of the quantum matter field through 
the Einstein-Langevin equation: 

G ab [g] = K((f ab \g}) + UMW (5.2) 

The stochastic source £, ab [g] at the- level of Gaussian approximation is completely characterized 
by its correlation function in terms of the noise kernel M abc d{x , y) ■ The noise kernel is the 
symmetrized connected part of the two-point quantum correlation function of the stress tensor 
operator with respect to the state of the matter fields. It describes their stress-energy fluctuations 
via 

{€ab\9> x )€cd\ff\y))t =N abcd (x,y) = - ({t ab [g;x),t cd [g]y)j} , (5.3) 

where t ab = T ab — (T a b) an d (. • .) is the usual expectation value with respect to the quantum 
state of the matter fields, whereas (. . .)g denotes taking the average with respect to all possible 
realizations of the stochastic source ^ ab . 

5.2.1. Open system approach: gravity interacting with quantum matter fields In the open 
quantum system approach one begins with a closed system, described by a density matrix 
p(x,q;x',q',t), defines a system S of special interest described, say, by the x variables or fields 
and look for the overall coarse-grained effects of its environments E, described by the q fields, on 
this system. The influence functional (IF) formalism |45] offers a very elegant way of capturing 
these effects on a microphysics level, with full self-consistency. If we are only interested in how 
the state of the system is influenced by the overall effects, but not the precise states, of the 
environment (s), then the reduced density matrix p r (x, x', t) = J dq p(x, q; x', q, t) could provide 
the relevant information. (The subscript r stands for reduced.) It is propagated in time from ti 
by the propagator J r : 



+00 +oo 

p r (x,x',t)= J dxi J dx'i J r (x,x',t | Xi,Xi,ti) p r (x i ,x' i ,t i ) (5.4) 

— oo — oo 

Assuming that the action of the coupled system decomposes as S = S s [x] + S e [q] + 
Sint[x,q], and that the initial density matrix factorizes (i.e., takes the tensor product form), 
p(x,q;x',q',ti) = p s (x,x' ,U)p e (q,q' the propagator for the reduced density matrix is given 
by 

J r (x,x',t\ x i ,x' i ,t i )= Dx / f Dx' ^ Bff \x,x'A 

J Xi J %'. 

where 

S e ff[x,x',t] = S s [x] - S s [x'] + S IF [x,x',t] (5.5) 

is the full effective action and Sif is the influence action. The influence functional T is defined 
as El 

T[x,x',t) = e lSlF[x ' x '' t] 

= [ dqfdqidqi F* Dq H Dq' ^ B [q]+S in *[x,q]-Seffl-S iut lz' pe ( q . } ^ t .y ( 5 g) 

' Jqi Jq'i 

Sif is typically complex; its real part TZ, containing the dissipation kernel D, contributes to the 
renormalization of S s , and yields the dissipative terms in the effective equations of motion. The 
imaginary part I, containing the noise kernel N, provides the information about the fluctuations 
induced on the system through its coupling to the environment. 

For a system with a quadratic action such as a simple harmonic oscillator, bilinearly coupled 
to a Gaussian environment (made up of e.g., non-interacting harmonic oscillators) as in the 
quantum Brownian motion (QBM) problem treated profusely in the literature (see, e.g., [461147]) 
the influence action takes the generic form : 

Sif [x, x'] = Z ■ A + A • H • E + • N • A, (5.7) 

where A = x' — x, E = |(x + x f ) and • = j/ dt. Z is a current term marking the contribution of 
the initial conditions, and the kernels H, N are related to the dissipation and noise kernels. In 
exact analogy with the above we can derive the dynamics of the gravitational field as an open 
system under the influence of N quantum matter fields viewed as its environment by calculating 
the influence action Sif defined as: (details are in [67], Appendix C) 

N r 

e iS 1F [h,h'] = pj / D0 i D^. e i5 »^'^- ift »^' fc ' ] /»[^ <) J <!>'/% (5.8) 



Up to quadratic order in the metric perturbations h a b it is given by 

Sif [Eo6, A afe ] = N ( Z ■ A + A • (H + M) • S + ^A • N ■ A ) , (5.9) 



5 For more details and further discussions on how the IF is applied in statistical field theory and applications, 
see, e.g., [66]. 



where we have introduced the semisum and difference variables T, a b = (h a j) + h' ab ) [2 and 
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(5.10) 
(5.11) 



where the notation T* was employed to indicate that the spacetime partial derivatives appearing 
in the time-ordered operators also act on the theta function implementing the time ordering. 
The functional derivative appearing on the right-hand side of (|5.1ip should be understood to 
account only for the explicit dependence on the metric: the implicit dependence through the 
field operator cj)[g] is excluded. 

As is known the kernels H and M exhibit divergences that are cancelled by renormalizing the 
gravitational coupling constant and the cosmological constant in the bare gravitational action 
as well as the coupling constants of the counterterms quadratic in the curvature. (We will not 
need terms of higher order in the metric perturbations because they give contributions to the 
connected part of the CTP generating functional of higher order in 1/N.) 



5.2.2. Closed system: semiclassical and stochastic gravity from large N expansion We now 
give a closed system treatment for the effects of an N component quantum matter field (f)j in 
a background spacetime with metric g ab . The semiclassical Einstein and the Einstein-Langevin 
equations can be identified as originating from the leading order and the next-to-leading order 
large N expansions. Treating perturbative quantum gravity (in terms of quantized linear metric 
perturbations h a b) interacting with N conformally invariant scalar fields (fij, Hartle and Horowitz 
|19] computed the effective action to leading order in large N and showed that the semiclassical 
Einstein equation can be interpreted as the equation governing the evolution of the expectation 
value of the metric to leading order [which for this theory is actually 0(1)]. The sources of the 
SCE equation are given by the expectation value of the stress tensor operator of the quantum 
matter fields. In the same vein Roura and Verdaguer [20] demonstrated that the next to leading 
order [which in this theory is 0(1/N)] [f| contribution to the quantum correlation functions in a 
large N expansion is equivalent to the stochastic correlation functions obtained in the context 
of stochastic semiclassical gravity. Let us see how this comes about. 

Consider metric perturbations around a globally hyperbolic background spacetime with 
metric g a b and scalar curvature R interacting with N minimally coupled free scalar fields 
= 1, ..,N). The action for the combined system is the sum of the gravitational action S g 
plus the action for the matter fields S m . The gravitational action is given by the usual Einstein- 
Hilbert term, the corresponding boundary term (which should be included to have a well-defined 
variational problem) and the usual counterterms required to renormalize the divergences arising 

6 Note the ordering scheme used in [201 167| is different from the customary usage. Their ieading order in the 
large N limit refers to the lowest order in 1/N with a nonvanishing contribution. Hence, as we will see, the 
leading order for the source of the semiclassical Einstein equation, which is proportional to the expectation value 
of the stress tensor operator, is 1/N°, or O(l) - this is called leading order large N in [19], whereas the leading 
order for the quantum two-point correlation functions is 1/N. In the counting scheme of [19] this would be the 
next-to-leading order which yields stochastic gravity. In the way how it is presented here this latter expansion is 
to be understood in a NLO sense. 



when functionally integrating the matter fields: 
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S=— d 4 xy^R(g) + — d^^K^ff) + (counterterms), (5.12) 

2k Jm K JS=dM 

where g a b = g a b + h a b is the perturbed metric, g a b is the background metric and the gravitational 
coupling constant k = %ir/m^ was rescaled to k/N so that the product of the rescaled 
gravitational constant times the number of fields remains constant in the limit N — > oo. The 
action for the matter fields is 

N -. 

S m = ~J2 d 4 x (r 6 V a ^V^. + m 2 tf) , (5.13) 

j=i jM z 

where m is the mass of the scalar field. We will use the expansion in 1/N to systemize the 
contributions. At the end one can always substitute back the rescaled gravitational constant in 
terms of the physical one. 

To make connection of the large N expansion with the theory of stochastic gravity one can 
use the procedure of Calzetta, Roura and Verdaguer [69J for the QBM problem to derive the 
symmetrized quantum correlation function for the metric perturbations (see Appendix C of |69|): 

ab 1 (i) 

-2 

+ ^ -N- (G TCt ) T ) abcd (x,x'). (5.14) 

There are two separate contributions to the two-point correlation function: the first one, called 
intrinsic fluctuations, is related to the dispersion of the initial state for the metric perturbations, 
whereas the second one, called induced fluctuations is proportional to the noise kernel and 
accounts for the fluctuations induced by their interaction with the quantum matter fields as 
its environment. Under the aforementioned conditions, the symmetrized quantum correlation 
function for the metric perturbations is equivalent to the stochastic correlation function obtained 
in stochastic semiclassical gravity by solving the Einstein-Langevin equation. 
The above expose serves to illustrate the following points: 

1) In a large N expansion semiclassical gravity (SCG) corresponds to a mean field theory 
obtainable from the leading order 0(1) expansion, and stochastic gravity (STG) theory 
corresponds to the next-to-leading 0(1/ N) order expansion. In this regard SCG and STG 
share the same properties as the 0(N) theory. Note however that stochastic gravity used in 
this more restrictive sense assumes Gaussian noise in the quantum matter field environment and 
weak perturbations off a background metric which is a solution of the SCE equations. In the 
more general sense such as used in [70] there is no such restriction and the parallel with the 
0(N) model ends at the Gaussian level. 

2) Stochastic gravity includes the fluctuations of the stress tensor in the matter sector via the 
noise kernel. This in turn enters into the determination of the two-point quantum correlation 
function for the metric perturbations. Contributions of higher order in 1/N enter either 
by including the vertices for the metric perturbations or terms from the influence functional 
evaluated beyond the Gaussian approximation, (e.g., in figure 1 of [67] the two diagrams shown 
give contributions of order 1/A^ 2 to the two-point quantum correlation function for the metric 
perturbations.) 

3) From these theories one can also see the relation between stochastic correlations and 
quantum fluctuations (via the loop expansion order), as well as the relation between correlation 
order in the Schwinger-Dyson hierarchy and the large N order. We shall dwell on the correlation 
aspects of MQP with the use of 2PI effective action in our sequel paper. 
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Note that it is by virtue of its Gaussian property that SCG as a mean field theory can be 
viewed as semiclassical. For STG if one does not make the Gaussian assumption for the noise, 
the Feynman- Vernon identity which enables one to interpret the quantum fluctuations of the 
matter field (as environment) as a classical stochastic source no longer holds. In what follows 
we give an example of field theories which are non-Gaussian. We will see that even at the 
rudimentary level the identification of mean field theory as classical is no longer valid. 

5.3. Large N gauge field and string theory 

5.3.1. 1/N expansion in gauge theories 't Hooft [24J introduced the 1/N expansion for 
nonabeliean gauge theories in order to express physical quantities as a systematic expression 
in powers of 1/N. Consider the gauge group SU(N). The gauge vector fields are Hermitian 
matrices A 1 ^, where both indices i,j run from 1 to N. The field strength is 

V, = dpAjj - d u A; 3 + ig[A^ A v % (5.15) 

where g is the gauge coupling constant. The Lagrangian density for a spinor particle of mass 
with wave function ip interacting with this gauge field is given by 

C = -~Tr(F^F^)-${i 7 D + m^i> 

= -Tr(i(«9 M A,) 2 - Id^dvAn + igd^A v {A^A u ] - ~g 2 [A„,A v ] 2 ) 

-${yfdp + ipfAp + m^. (5.16) 

It is convenient to introduce a double line notation to keep track of the group indices i and 
j. Whenever there is a close index loop in the Feynman diagram, we get a factor of N from 
counting the values of the index can assume. This is how the dimension N of the gauge group 
enters. Since all three-point vertices come with a factor of g, and all four-point vertices come 
with a factor of g 2 , the g and N dependence of the amplitude for a general Feynman diagram is 
seen to be (g 2 N) F+2x ~ 2 N 2 ~ 2x , where F is the number of loops and x is the genus, or number 
of "holes" of the polyhedron formed by the Feynman diagram. In the 't Hooft large N limit 
while keeping g 2 N = A fixed, we see that the N dependence of the diagram is N 2 ~ 2x . So the 
diagram without holes (x = 0) dominate. In other words, in this limit, only the planar diagrams 
contribute. The interactions between hadrons is expected to be a 0(1/N) effect. Hence in the 
large N limit one can consider the problem of confinement and hadron mass spectrum without 
the complication of residual hadron interactions. The hope was that one could solve the theory 
with N = oo exactly followed by an expansion in 1/N = 1/3 for QCD. 

Witten [25] referred to this limit as analogous to the classical limit for quantum mechanics, 
in the sense that for large N there exists one single "classical" gauge field, called "master 
field", which saturates the path integral and gives the dominant contribution. Such a classical 
field configuration was called a Witten's argument was based on the assumption that all possible 
invariant operators satisfy the factorization property [which holds in any U(N) symmetric theory 
of N x N matrix fields Aij(x)} 

k 

(/ 1 (A).../ fe (A)) = n(/i(^))[l + C(l/iV)], (5.17) 

8=1 

where the fi(A) are U(N) symmetric operators of the form f(A) = Tr[A(xi)...A(xi)]. 

However, Haan [32J pointed out that there is a gap in Witten's argument for the existence of 
a master field in the large N limit of QCD. He showed that the master field does not exist for 
the two-matrix model and the large N limit is not a classical limit. Instead, the factorization of 
expectation values of invariant operators is shown to be a property analogous to the vanishing 
of fluctuations for macroscopic observables in the thermodynamic limit. 



5.3.2. Large N field theories and string theory At large N, QCD is expected to behave like a 
string theory. This is supported by the following observations. First, Feynman diagrams are 
organized in a topology expansion, just like string worldsheets. Second, the existence of Regge 
trajectories indicates a string-like behaviour. Third, QCD contains string-like objects, which are 
the electric flux tubes between quarks and antiquarks. Maldacena's AdS/CFT correspondence 
|26j conjectures the equivalence of A = 4 SU(N) super Yang-Mills theory with gauge coupling 
gYM in 4 dimensions and Type IIB superstring theory with string coupling g s in AdS^ x S* 5 
space where both AdS§ and 5 5 have the same radius L. The parameters between these two 
theories are related by 

g s = g YM , L 4 = A7T 9s N(a') 2 , (5.18) 

where a' is related to the string tension. This conjecture is assumed to be valid for all values of 
N and of g s = g\ M . Assuming the validity of Maldacena's conjecture, one can study different 
limits of field theory and string theory. 

• First take the 't Hooft limit. This is the case where A = 5y M A = g s N fixed and letting 
N — > oo. In the Yang-Mills theory, this limit is well-defined, at least perturbatively, and 
corresponds to a topological expansion of the field theory's Feynman diagrams. On the 
AdS side, the string coupling is g s = X/N and the 't Hooft limit corresponds to weak 
coupling string perturbation theory. The 1/N expansion in gauge theory corresponds to 
the g s string loop expansion in the string theory side. Hence the Maldacena conjecture in 
this case becomes a correspondence between large N limit of gauge theories and the classical 
string theory. (Here classical means there is no string loop expansion in the N — > oo limit). 

• The large A limit. In this limit, one can expand the effective action in the AdS side in terms 
of ol . The distance scale is set by the AdS radius L which is related to the scale of Riemann 
tensor by R ~ 1/L 2 = (g s N)~2 jol = A _ 2 jol . Hence in the large A limit, the classical Type 
IIB string theory on AdS^ x S 5 becomes classical Type IIB supergravity on AdS§ x S 5 
and the a' expansion in string theory corresponds to A~2 expansion in field theory. (Here 
gravity is weak because the curvature scale approaches zero in this limit and one can neglect 
the stringy correction and approximate it by the classical supergravity action). 

To summarize, let us see what all this means for the issues we wish to address in this paper: 
classicality and large N from the perspective of string/gauge/gravity theories. 

Old school field theorists often refer to the large N limit of a Yang-Mills gauge theory as 
classical. Or stretching it a little, one may habitually view all large N theories as classical. The 
advent of string theory alters fundamentally the way we look at nature, not just pertaining to 
Planck scale microphysics but also in relation to gauge fields and gravity, two essential elements 
of nature as we see it today. To begin with, N enters in a new coupling constant called the 
't Hooft coupling, in addition to the string tension, the magnitude of both and their interplay 
alter a theory's qualitative behaviour in a fascinating way. Another important new concept 
introduced by string theory is duality: between string theory and gauge theory, and between 
gauge theory and gravity. 

The correspondence between gauge theory and string theory provides another window to 
examine the large N behavior of gauge theory. Under this duality relation, a Yang-Mills gauge 
theory ingrains also the effects of gravitational interaction: Under different regimes of 't Hooft 
coupling, the behavior of large N limit seen from the AdS side of the correspondence can 
appear either as the suppression of string loop expansion and hence a classical string theory or 
the suppression of the stringy effect giving rise to a weakly interacting gravity theory where the 
quantum effect is negligible. 



We can only give an inkling of the immense body of work accumulated in the last half a 
century for large N, especially since 1998 when the AdS/CFT correspondence was proposed. 
The simple message we wish to convey here is: "exercise caution when tackling the meaning or 
finding the criteria of classicality" : the explanation of the emergence of classical behavior or even 
the apparently straightforward identification of a classical regime from an interacting quantum 
theory is a subtle issue, as it may contain multiple structures and carry different dynamics in 
different interaction regimes. 

6. Conclusion 

Macroscopic quantum phenomena is a relatively new research venue, with exciting experiments 
and stimulating proposals. Surprisingly there are very little systematic theoretical investigations 
or even serious attempts we are aware of to construct a viable theoretical framework for this 
new endeavour. In this paper we ask only one simple question: Is it really that surprising 
that quantum features can appear in macroscopic objects? Note the surprise comes from the 
conventional and hitherto unchallenged view that macroscopic means classical. The purpose 
of our series of papers is to remove this degeneracy. By examining many well-known examples 
where theoretical analysis has been performed in detail for different systems, we find that there 
is no a priori good reason why quantum phenomena in macroscopic objects cannot exist. The 
challenge for theorists is to identify for various systems the physical conditions MQP manifest 
and in what ways. 

In this paper we take the large N perspective, and try to understand what the leading order 
and next-to- leading order large N expansions bring us. In general we can only say that the 
LOLN gives mean field theory and the NLO gives deviations from the mean field, both are 
still very much quantum in nature. Only for Gaussian models such as the 0(N) theory (zero- 
dimensional field theory, referring to the quantum mechanics of N oscillators) or in theories where 
the Gaussian approximation is imposed, such as semiclassical gravity and stochastic gravity, do 
we find that the mean field at LOLN is equivalent to classical, and the NLO incorporates 
deviations from the mean, which, only when one identifies the mean with classical would it 
be conjured as containing "quantum fluctuations". By combing through these examples we 
are trying to change the common attitude on a very subtle point, namely, "Don't begin with 
classical and try to go backwards to quantum. One should always begin with quantum and 
only under very special conditions and stringent criteria would a quantum system begin to show 
certain classical features." We want to shift the burden of proof to those who believe they live 
in a classical world, to find these conditions and check every criteria. If they are not completely 
fulfilled, well then, perhaps one should just accept that the world is fundamentally quantum 
and it is only under specific yet commonly encountered conditions H that it appears classical. 

This message is of course not that foreign to researchers working in the last two decades on 
the quantum to classical transition or the Q-C correspondence issues. The half a dozen or so 
criteria we reviewed in this paper stem from these investigations. What we have done somewhat 
different here is to examine every criteria anew with a macroscopic system, here with the aid 
of large N expansions. There is an abundance of information in the particle physics literature 
on this subject. We brought forth some here to serve as illustrations for the one question we 

7 Maybe this is a sociological phenomenon, that exciting experiments are carried out by atomic-optical physicists 
while those who care or know more about the theories are mainly in the particle and gravitational physics 
community. If so this series of papers can be viewed as a modest attempt to bridge that divide. 

8 Specific and common are not contradictory conditions: equilibrium state is an example. It is one specific state 
out of a large number of states accessible to a large system, but it is also the most commonly encountered - with 
high probability - when the system is placed in a heat bath. Behind this is the working of central limit theorem, 
quantum regression theorem in counting the statistics, or the conservation laws governing the hydrodynmic 
variable, as far as decoherence is concerned. These theorems and laws give a deeper meaning of 'habitual'. 



raised, but there is plenty more material which can be used fruitfully for performing quantitative 
analysis of MQP. 

There are other important issues which we have left out. For example, here we use the number 
of components N as a measure of macroscopic, but it could be physically more meaningful to 
incorporate interaction strengths and dynamics, in the renormalization group sense. Thus it 
could be the ultraviolet versus the infrared behaviour of a (quantum) system which conveys a 
stronger microscopic versus macroscopic sense, as manifests in critical phenomena. It is near 
the critical point that the system's IR behaviour dominates, where the correlation function 
extends to infinity. This aspect of correlation, both classical and quantum, is the focus of our 
sequel paper. For interacting quantum systems we will invoke the nPI effective action and the 
Schwinger-Dyson equations (corresponding classically to the BBGKY correlation hierarchy) to 
analyze what quantum versus classical means for large/small and/or strongly/weakly correlated 
systems. In that context we will include critical phenomena and other statistical mechanics 
considerations to examine this new issue of MQP at the foundation of quantum mechanics. 
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